then 
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Two easily applied inequalities are given for the 'Hail probabilities" of multivariate 
normal distributions. 

A basic integral arising in statistics is of the form: 

i]y[|i/2 r^ r~ r 1 "1 

i^(a,M)=|gL_J^ ...J^ exp|^-i(x+a)M(x+a)'Jr/.ri...r/x, (1) 

where a and x are n component row vectors, M is an nXn positive definite symmetric matrix 
with determinant |M| and '^'^' n.ieans transpose. In particnhir (1) arises in evahiating 
Pr(X>b) where X is a random variahk* wit li a nudtivariate normal distribution, and b is a vector 
of constants. In this note iiicfjuahties (up])er and lowei- })Ounds) are^Jpresented for F{a,M). 
Let 

7<^(a,M)=/(a,M) r ... f exp (-aMxOXexp ("-ixMx'y/ri . . . r/r,. (2) 

Assimie throughout tlie following, 

A=aMyO/iA\, Ai=^ ajrrijiyo iov i = \ , . . .,n. (3) 

The following elenientar}' inequalities are required 

2r+l 2r 

J2 (-cy/k\<:e-'<: J2 {-cy/k\ for r=0, 1 , . . . (4) 

and c>0. In particular, (4) will l)e used with c = )^xMx', the last quantity being positive 
whenever xf^O, i.e., at least one Xi^^O, since M is positive definite. 
Now using the right-hand side of (4) with r=0 in (2) one obtains 

7^(a,M)</(a,M) J^" . . ,j^^ exp (~-|] A,r,) dx„ . . ., ./.r„. (5) 

After evaluating the integral in (5) one has 

i^(a,M)</(a,M)(_nA,)"'. (I) 

Now using tlie lel'l-liaiul side of (4) with r=0 in (2) one obtains 

F(a,M)>/(a,M) |*^ . . . J^ (l-^f] S rinnix^X exp (^-g A,^,^ dx„ . . ., dx„. (6) 
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After evaluating the integral in (6) one obtains 

F(a,M)>/(a,M) ( n A,)"' [l-g S 2 m,/l+5,,)/(A4,)]. (II) 

where dij=l if i=j and 8ij=0 if tT^j- 

Inequalities similar to (I) and (II) can be obtained by using r>0 in (4). These inequalities 
would be complicated notationally. In particular, however, Feller [1, pp. 166, 179]^ presents 
these inequalities (for general r) when n=l and mu = l. Polya [2] has considered the case of 
r>0 for the integrand in (1) for a different region of integration. Polya's sequence of inequali- 
ties has the desirable property of improving as r increases until the sharpest form is obtained 
and then of becoming successively weaker. Ruben [3] has considered the special case where 
mii = l for i=l, . . ., n and mij=m for i^j. His results appear to be complicated. 

If the case considered by Ruben is specialized in the following manner: 7i=2 and ai=a2=a, 
then one obtauis from using (4) in (2) the following 

Fia,M)m^,M) Ij AH-2A^)- [s'" ^'"+^'^'1,:;+^^'^'"^° ] (HI) 

where > (<) if i* is odd (even), A=l + m, and 2^'^ is over all u,v,w>0 such that u + v+iv = i. 
It is not clear that the inequalities in (III) have the desirable feature of the Polya bounds. 
Example 1. Assume X has a bivariate normal distribution with mean vector zero, and 



variance co variance matrix M ^- 



1 1/2 
1/2 1 
pression for the desired probability is 



Approximate Pr(X>(3,3)). The integral ex- 



j I f{x,M)dxidx2 



In this integral make the following change of variables yi=Xi—?> and y2=X2—3. The integral 

4/3 -2/3 
2/3 4/3 



then becomes F(a,M) where a==(3,3), M= 
Now using (I) one obtains 



and !M|=:4/3. Note Ai = A2=-2. 



Pr(X>(3,3))< W-3; =lq.qoqii4^ 

ox 

and from (II) one obtains 

From [5] one finds the correct value to be 0.000082. 

In this example the ratio of the upper to the lower bound is 2. If one considers the same 
problem, with positive correlations of 1/2 replaced by —1/2 the corresponding ratio is 54/49. 
Since the probabilities are smaller with negative than with positive correlations, this is not a 
surprising result. 



4/3 2/3 
2/3 4/3 



Example 2, Assume M= 

then 

0.0000039>i<^(a,M)>0.0000031 



so that |M|-=4/3, also assume a= (2,2) so that A= (4,4) 



5 Figures in brackets indicate the literature references at the end of this paper. 
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To (*()m|)ar(' tlu* acciu'iicN' of the results lor (ho biviiriate case witli the univariate case, consider 
the foUovvin^- univariate example where the true value is about the same as in the above bivariate 
examj)l(\ Let 7i=^l, 7nn~^, and a = 4.5, then one obtains 

0.0000()355>/^X4.5, 1)>0.00000337. 

The correct value, from [4], is 0.00000339S. 

In the present notation, the genei'al result foi' the univariate case given by Feller [1, p. 179] 
is 

F{a,\)^a-^f{a,\)j: ^ W ' 

where a>0 and >(<) if R is odd (even). 

Exainple 3. If ma=l for ^=1, . . ., n and 7njj = m for i^j, and Si=(a, . . . , a), then |M| = 
(l—my~'^[l+(n—l)m] and Aj=a[l + (n— l)//^]. (To preserve the positive definite character 
of M it is necessary that l + (7i— l)m>0, i.e., m> — (n— 1)"^ and 7n<^l.) Now using (1) and 
(IT), one obtains 

/^Xa,M)<|M|'/^f^-"'^^/-(27r)-"/^A-'^ 
and 

/^(a,M)>|M|^/^V-^''^^/-(27r)-"/'A-"ri-^'^!^+'^--^^^ 

where A=l + (/i— l)m. ( 'onsider the case with it=l], m=^l/3 and a = o{l]/2)^^-. (Tiiis numeri- 
cal example corresponds to finding the probability that a tiivariate normal variable with mean 
vector zero, and covariance matrix having I's on the diagonal and 1/2's off' the diagonal has each 
component larger than 5.) In this case the ratio of the upj)er })ound to the lower bound 
of 7^(a,M) is 25/13. If one considers a = 9{3/2y^'^ the corresponding ratio is 27/23. 

With the previous notation and restrictions, considei- a new integral, 6^(a,M) 
obtained from modifying F{si,Wi) by replacing the irpper limits of integration by ones, i.e., 

G{a,M)= r . • r'/'(x+a,M)r/.r, . . . r/./v (7) 

Jo Jo 

Then evaluating the integral analogous to (0), one obtains , 

«(a,M)</(a,M) ( n|l-<'-^'l/A,)- (IV) 

From the analogue of (0), one obtains 
a(,,M)>/(a,M)[n (!-,-.,/.,] {, -1 [S ^ (l + '-';^!i7'" ) 

Results (IV) and (V) w^ere obtained by Professor Olkin. 



A first draft of this paper was prepared in 1954 while the author was at the National Bureau 
of Standards. That work was stimulated by Dr. Benjamin Epstein in conversations on the 
distribution of extreme values for dependent variables. Professor Ingram Olkin suggested 
changes in the original draft which have been incorporated at this time. During the sunnner 
of 1961, under a National Science Foundation Undergraduate Research Program, Thomas 
Kieran did some additional numerical analysis for the manuscript. 
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